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The Boltzmann equation for inelastic Maxwell models is used to determine the Navier-Stokes 
transport coefficients of a granular binary mixture in d dimensions. The Chapman-Enskog method 
is applied to solve the Boltzmann equation for states near the (local) homogeneous cooling state. 
The mass, heat, and momentum fluxes are obtained to first order in the spatial gradients of the 
hydrodynamic fields, and the corresponding transport coefficients are identified. There are seven 
relevant transport coefficients: the mutual diffusion, the pressure diffusion, the thermal diffusion, the 
shear viscosity, the Dufour coefficient, the pressure energy coefficient, and the thermal conductivity. 
All these coefficients are exactly obtained in terms of the coefficients of restitution and the ratios of 
mass, concentration, and particle sizes. The results are compared with known transport coefficients 
of inelastic hard spheres obtained analytically in the leading Sonine approximation and by means 
of Monte Carlo simulations. The comparison shows a reasonably good agreement between both 
interaction models for not too strong dissipation, especially in the case of the transport coefficients 
associated with the mass flux. 
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I. INTRODUCTION 

The evaluation of the transport coefficients from the Boltzmann equation for inelastic hard spheres (IHS) is quite 
involved. In fact, to get explicit results one usually considers the leading order in a Sonine polynomial expansion of 
the velocity distribution function. These difficulties increase when one considers multicomponent systems since not 
only the number of transport coefficients is larger than for a single gas but they are also functions of more parameters 
such as composition, masses, sizes, and the coefficients of restitution. As in the elastic case, a possible way to partially 
overcome these problems is to consider other interaction models that simplify the complex mathematical structure 
of the Boltzmann collision integrals for IHS. For this reason, the so-called inelastic Maxwell models (IMM) have 
been widely used in the past few years as a toy model to characterize the influence of the inelasticity of collisions 
on the physical properties of granular fluids. The IMM share with elastic Maxwell molecules the property that 
the collision rate is velocity independent but, on the other hand, their collision rules are the same as for IHS. In 
this sense, although these IMM do not correspond to any microscopic potential interaction, it has been shown by 
several authors 1 ' 2 ' 3 ' 4 ' 5,6,7,8 ' 9 that the cost of sacrificing physical realism is in part compensated by the amount of exact 
analytical results. 

Most of the studies carried out by considering IMM have been devoted to homogeneous states, especially in the 
analysis of the overpopulated high energy tailsA&LSilfiiAi However, much less is known in the case of inhomogeneous 
situations and, more specifically, on the dependence of the transport coefficients on dissipation. For a monocomponcnt 
granular gas subjected to simple shear flow, the IMM have been used to calculate the rheological properties (shear and 
normal stresses) in three dimensions; 12 More recently, this study has been extended 1 ^ to multicomponent systems, 
the exact results of IMM showing a close agreement with those obtained analytically for IHS in the first Sonine 
approximation 14 and by means of Monte Carlo simulationsi^ii 5 . All these results are restricted to steady shear flow 
problems without any limitation on the strength of the shear rate. For general inhomogeneous problems and in the case 
of a monocomponent gas, the Boltzmann equation for IMM has been solved 1 ^ from the Chapman-Enskog methodic for 
states near the (local) homogeneous cooling state. Explicit expressions of the Navier-Stokes transport coefficients of 
IMM in d dimensions have been obtained for unforced systems as well as for systems driven by thermostats. In contrast 
to the findings of Ref. Il3l the comparison with the transport coefficients of IHS 18,19 shows that their dependence on 
inelasticity is captured by the IMM only in a mild qualitative way. This fact stimulates the determination of the exact 
expressions of the transport coefficients for inelastic granular mixtures. 

The goal of this paper is to derive the hydrodynamic equations for a d-dimensional binary mixture of inelastic 
Maxwell gases at low-density. As in the single gas case^ a normal solution to the coupled set of Boltzmann equations 
for the two species is obtained by using the standard Chapman-Enskog methodic conveniently generalized to inelastic 
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collisions. In the first order of the spatial gradients of the hydrodynamic fields, we get the corresponding Navier-Stokes 
hydrodynamic equations with explicit expressions for the relevant transport coefficients of the mixture. For molecular 
fluid mixtures, the specific set of gradients contributing to each flux is restricted by fluid symmetry, time reversal 
invariance (Onsager relations), and the form of entropy production^ However, in the case of inelastic collisions, only 
fluid symmetry holds and so there is more flexibility in representing the fluxes and identifying the corresponding 
transport coefficients. In the case of the pressure tensor, fluid symmetry considerations implies that its form to first 
order in the gradients is the same as for the monocomponent gas. In the case of mass and heat fluxes, several different 
(but equivalent) choices of hydrodynamic fields can be used and some care must be taken when comparing transport 
coefficients coming from different representations. As in the case of IHS^Si we take the gradients of the mole fraction, 
the (hydrostatic) pressure, the temperature, and the flow velocity as the relevant ones. As a consequence, in this 
representation there are seven independent scalar transport coefficients: the mutual diffusion, the pressure diffusion 
and the thermal diffusion associated with the mass flux, the shear viscosity corresponding to the pressure tensor and 
the Dufour coefficient, the thermal conductivity, and the pressure energy coefficient associated with the heat flux. All 
these coefficients are given in terms of the coefficients of restitution as well as on the ratios of concentration, masses and 
particle sizes. In addition, as in the previous study for IHSj 2 ^ our theory takes into account the effect of temperature 
differences (failure of energy equipartition) on the transport coefficients, leading to additional dependencies of them 
on the concentration. 

The plan of the paper is as follows. In SectionlTTIthe Boltzmann equation for IMM and the macroscopic conservation 
laws are introduced. The model includes average collision frequencies uo rs which can be freely fitted to get good 
agreement with IHS. Here, we fix uj rs by the criterion that the cooling rates £ rs of IMM be the same as those 
obtained for IHS in the local equilibrium approximation. The homogenous solution of the Boltzmann equation is 
analyzed in Section IlIII where the temperature ratio and the fourth cumulant (kurtosis) of the velocity distribution 
functions are exactly obtained. Comparison with the results obtained for IHS 22,23 shows an excellent agreement 
for the temperature ratio but significant discrepancies with the fourth cumulant of IHS. Section llVI deals with the 
application of the Chapman-Enskog method to get the transport coefficients of IMM. In Sectional the dependence of 
some of these coefficients on the parameters of the system is illustrated and compared with known results derived for 
IHS. 22 24-25,26 The comparison shows in general a qualitative good agreement, especially for the transport coefficients 
defining the mass flux. The paper ends in Section IVII with a brief discussion on the results reported in this paper. 



II. INELASTIC MAXWELL MODELS FOR A GRANULAR BINARY MIXTURE 

Let us consider a binary mixture of inelastic Maxwell gases at low density. The Boltzmann equation for IMM 3 ^* 1 ^ 2 ^ 
can be obtained from the Boltzmann equation for IHS by replacing the rate for collisions between particles of species 
r and s by an average velocity-independent collision rate, which is proportional to the square root of the "granular" 
temperature T. This means that a random pair of colliding particles undergoes inelastic collisions with a random 
impact direction. With this simplification, the velocity distribution functions f r (r,v;t) (r = 1,2) of each species 
satisfy the following set of nonlinear Boltzmann kinetic equations: 

(0 t + v-V)/ r (r,v;i) =5Z J «M/r(*). /.(*)] > (1) 

s 

where the Boltzmann collision operator J rs [v|/ r ,/ a ] is 

J rs vi / r ,/ s = — / dv 2 / da 

n s [r,t)il d J J 

x [a i : s 1 / r (r,v' 1 ,i)/ s (r,v^,t)-/ I .(r,v 1 ,i)/ s (r,v 2 ,t)] . (2) 

Here n r is the number density of species r, uj rs ^ co sr is an effective collision frequency (to be chosen later) for 
collisions of type r-s, Qd — 2n d / 2 /T(d/2) is the total solid angle in d dimensions, and a rs = a sr < 1 refers to the 
constant coefficient of restitution for collisions between particles of species r with s. In addition, the primes on the 
velocities denote the initial values {v' 1; v 2 } that lead to {vi,v 2 } following a binary collision: 

vi = vi - fi sr (l + a~}) (a ■ gi 2 )er, v 2 = v 2 + [i rs (l + ot~}) (a ■ g^)? , (3) 

where gi 2 = vi — v 2 is the relative velocity of the colliding pair, <x is a unit vector directed along the centers of the 
two colliding spheres, and \i rs = m r /(m r + m s ). The collision frequencies to rs can be seen as free parameters in the 
model. Its dependence on the restitution coefficients a rs can be chosen to optimize the agreement with the results 
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obtained from the Boltzmann equation for IHS. Of course, the choice is not unique and may depend on the property 
of interest. 

There is another more refined version of the inelastic Maxwell modeli*SiS£ where the collision rate has the same 
dependence on the scalar product (<r • §12) as in the case of hard spheres. However, both versions of IMM lead to 
similar results in problems as delicate as the high energy tails^ Therefore, for the sake of simplicity, here we will 
consider the version given by Eqs. (0 and ©. 

The relevant hydrodynamic fields in a binary mixture are the number densities n r , the flow velocity u, and the 
granular temperature T. They are defined in terms of the distribution f r as 

dv/ P (v), (4) 

PU = ^2 Pr u r = ^ J dvm r vf r (v), (5) 

r r 

nT = p = J2 n r T r = ]T ^ J dvV 2 f r (v), (6) 

r r 

where p r = m r n r is the mass density of species r, n = n\ + n 2 is the total number density, p = p\ + pi is the total 
mass density, V = v — u is the peculiar velocity, and p is the hydrostatic pressure. Furthermore, the third equality 
of Eq. (jSJ defines the kinetic temperatures T r of each species, which measure their mean kinetic energies. 

The collision operators conserve the particle number of each species and the total momentum, but the total energy 
is not conserved: 

J dvJ rs [v\f r ,f s ]=0, (7) 
/"rfvm r vJ rs [v|/ r ,/ s ] =0, (8) 

r,s J 

I dv^m r V 2 J rs [v\f r J s ] = - d -nTC, , (9) 

r,a J 

where £ is identified as the "cooling rate" due to inelastic collisions among all species. At a kinetic level, it is also 
convenient to introduce the "cooling rates" Q r for the partial temperatures T r . They are defined as 

Cr - ]T Crs = - ]T f dvm r V 2 J rs {v\f r , f s ], (10) 

s s r r 

where the second equality defines the quantities £ rs . The total cooling rate £ can be written in terms of the partial 
cooling rates Cr as 

C = T- 1 Y, X rT r Cn (11) 

r 

where x r — n r /n is the mole fraction of species r. 

From Eqs. (@J) to ©, the macroscopic balance equations for the binary mixture can be obtained. They are given 
by 

D t n r + n r V-u-\ — =0 , (12) 



DtU + p^'V • P = , 



(13) 
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D t T V — 

71 — ' m.~ 



dn 



(V • q + P : Vu) 



-CT 



In the above equations, D t = dt + u • V is the material derivative, 

jr = m r J dvV/ r (v) 
is the mass flux for species r relative to the local flow, 

P = J2 f dvTO r VV/ r (v) 

r 

is the total pressure tensor, and 

q = E/ dv±m r V 2 Vf r (v) 



(14) 



(15) 



(16) 



(17) 



is the total heat flux. The balance equations (|12|) - (|14[> apply regardless of the details of the model for inelastic 
collisions considered. However, the influence of the collision model appears through the dependence of the cooling 
rate and the hydrodynamic fluxes on the coefficients of restitution. 

As happens for elastic collisions^ the main advantage of using IMM is that a velocity moment of order k of the 
Boltzmann collision operator only involves moments of order less than or equal to k?^ This allows one to determine 
the Boltzmann collisional moments without the explicit knowledge of the velocity distribution function. The first few 
moments of the Boltzmann collision operator J rs [f r , f s ] have been explicitly evaluated in Appendix IA1 In particular, 
according to Eq. 1)1(1 fl . the second moment of J rs [fr, fs] allows one to get a relationship between the collision frequencies 
uj rs and the cooling rates Crs- From Eq. I|A2(1 . one easily gets 



2uj r 



-fi sr (l + a rs ) 



1 



■(1 + a rs )- 



H sr (l + a rs ) - 1 . . 

+ j Jr ' Jf 



where 



(18) 
(19) 



p r = n r T r is the partial pressure of species r and j r = T r /T . 

In order to get explicit results, one still needs to fix the parameters lo ts . The most natural choice to optimize the 
agreement with the IHS results is to adjust the cooling rates ( rs for IMM, Eq. (|18fl . to be the same as the ones found 
for IHSi2£ Although the cooling rates are not exactly known for IHS, a good estimate of them can be obtained by 
considering the local equilibrium approximation for the velocity distribution functions / r , i.e., 



/r(V) 



2-KT r 



d/2 



exp 



In this approximation, one has22 

*IHS 



+ 1 



1/2 



(1 + a rs ) 



m r V z 
2T r 



-. rsr \ ■ 

1-— (l + a rs )- 



(20) 



(21) 



where vo{t) = \/2T{m\ + mi)jm\mi is a thermal velocity defined in terms of the temperature T(t) of the mixture. 
Thus, according to Eq. (|18fl . the collision frequencies u> rs are given by 



co rs = Ax. 

where vq is an effective collision frequency given by 



CT12 



d-1 



1/2 



40F ? 



4~w 



(22) 



(23) 



Upon deriving (|22|l use has been made of the fact that the mass flux j r vanishes in the local equilibrium approximation 
(|2Up. In the remainder of this paper, we will take the choice (1221 for Lo rs . 
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III. HOMOGENEOUS COOLING STATE 

As a previous step to determine the Navier-Stokes transport coefficients from the Chapman-Enskog method^ 
one needs to analyze the homogeneous solution of the Boltzmann equation l|T]l. In this case (spatially isotropic 
homogeneous states), the Boltzmann equation (JTJ becomes 

d t f r (v,t) = J2Jrs[v\frJs]. (24) 
s 

From Eq. 1(24(1 . one has 

d t T = -CT, d t T r = ~( r T r . (25) 
The time evolution of the temperature ratio 7 = T\(t) /Ta(i) follows from the second equality of Eq. (|25|) : 

9 t ln 7 = C 2 -Ci- (26) 
In addition, since the mass flux vanishes in the homogeneous state, Eq. ((18(1 gives the relation 



Cr — ^ Crs — ^ , f^sr (1 ~t~ C^rs) 



— (1 + <*„)- 



(27) 



where uj rs is given by Eq. (|22(l . 

The so-called homogeneous cooling state (HCS) qualifies as a normal solution for which all the time dependence of 
f r {v\t) is through the global temperature T(t). Consequently, it follows from dimensional analysis that f r (v;t) has 
the scaling form 

f r (v, t) = n r Vn d (t)<P r (v/v (t)), (28) 

where Vq is the thermal velocity previously introduced. The fact that the time dependence of f r (v,t) only occurs 
through the temperature T{t) (which is the relevant one at a hydrodynamic level) implies that the three temperatures 
Ti(t), T2{t) and T(t) are proportional to each other and their ratios are all constant. One possibility is that all three 
temperatures are equal, as in the case of elastic collisions. However, as we will see later, the partial temperatures are 
different. Since the temperature ratio must be independent of time, Eq. ((26(1 leads to the equality of the cooling rates: 



Ci(t) = C2(*)=C(*). (29) 

When the expression (|27|l is substituted into 1)29(1 , one gets a closed nonlinear equation for 7 whose numerical solution 
gives the dependence of the temperature ratio on the parameters of the problem. Except for mechanically equivalent 
particles, our results show that although both species have a common cooling rate, their partial temperatures are 
clearly different. This implies a breakdown of the energy equipartition. The violation of energy equipartition in 
multicomponent granular systems22i2Liffli2ii22i2ii has also been confirmed in computer simulations^^ii&Sii and even 
observed in real experiments in two^£ and three^ dimensions. It must be remarked that the fact that Ti(t) 7^ ?2(t) 
does not mean that there are additional hydrodynamic degrees of freedom since the partial temperatures T r can be 
expressed in terms of the granular temperature T as 

I + £1(7- I) 1 + 0:1(7-1) 

The problem is therefore to solve the Boltzmann equation for a distribution of the form l|28|l subject to the self- 
consistency constraint (|29|l . In terms of the reduced velocity v* = v/vq, the Boltzmann equation (|24|) for the reduced 
distribution $ r defined in Eq. 128(1 becomes 



5C^-(v*^)=E J "[«1^* S ], (31) 



with = C/^o- I n addition, 

J?.[*r, **] = ^ J dv* J da [«-/<&,(<)$,(<) - $ r K)$sK)] , (32) 




O 1 — ' 1 1 ■ ' ■ 1 1 1 

0.5 0.6 0.7 0.8 0.9 1.0 



« 

FIG. 1: Plot of the temperature ratio 7 = T1/T2 versus the coefficient of restitution a for d = 3, £1 = |, ai/02 = 1 and three 
different values of the mass ratio 7711/7712: (a) 7711/7712 = 10 (circles), (b) mi/7712 = 2 (squares), and (c) 7711/7712 = 0.1 (triangles). 
The solid lines are the results derived here for IMM, the dashed lines correspond to the results obtained for IHS from the first 
Sonine approximation, 22 and the symbols refer to Monte Carlo simulations for IHSj— 



where uj* = uj r /i>Q. Upon writing Eq. I|31|) . we have accounted for the time dependence of f r , which implies that 

d t f r = ~(Td T f r = US- ■ (v/ r ) . (33) 
2 ov 

Although the exact form of the distribution <I> r is not known, an indirect information on the behavior of $ r is given 
through its velocity moments. In particular, the deviation of $ r with respect to its Maxwellian form <I> r ,M can be 
characterized through the fourth cumulant 



where 



and 



r V 4 },-l 



d{d+2 r 

dv*V* 4 ^ r (v 



<f> rtM (v*)=ir- d / 2 6f /2 e- erV * 



(34) 



(35) 



(36) 



To get the fourth velocity moment, we multiply both sides of the Boltzmann equation 1)31(1 by v* A and integrate over 
the velocity. The result is 



(37) 



The collisional moment appearing on the right-hand side of l|37|l is evaluated in Appendix ^] with the result 

li sr (l + a rs ) 



J dv*«* 4 j;j$ r , $ s 



iM+ ^ s {3/^(1 + «, S )V 4 }, 

+ [2d + 3/4(1 + a rs f - 6/i sr (l + a rs ) + 4] [ Msr (l + a rs ) - 2] (v* 4 ) r 
+ d ( d + 2 ) /Jar g-ig-i( 1 + ars) [2d + 4 - 12/x sr (l + a rs ) 
+6fi 2 sr (l + a rs ) 2 ]} . 



(38) 



Substitution of Eq. I|38|) into Eq. (|37|l leads to a coupled set of linear equations for (f* 4 )i and (w* 4 )2 (or equivalently, 
for ci and C2). The solution of this set gives ci and ci in terms of the parameter space of the problem. In the 
one-dimensional case (d — 1), our results reduce to the ones previously obtained by Marconi and Puglisi 32 from the 
so-called Maxwell scalar model (i.e., by taking oj rs oc n s ). Moreover, for mechanically equivalent particles (cti = a-2 1 
mi = 7Ti2, an = (^22 = ^12 = a), the results obtained by Santos 1 ^ for the single gas case are recovered, namely, 7 = 1 
and 



ci = c 2 



12(1 -a) 5 



Ad- 7 + 3a(2 - a) 



(39) 



7 



0.5 - 




a 



FIG. 2: Plot of the coefficients c r versus the coefficient of restitution a for d = 3, xi = |, (J\ja% = I and m\jmi = 2. The solid 
lines are the results derived here for IMM, the dashed lines correspond to the results obtained for IHS from the first Sonine 
approximation;— and the symbols refer to Monte Carlo simulations for IHSi—The circles correspond to ci while the squares 
correspond to C2. 



A full presentation of the results is difficult due to the many parameters involved in the problem: an, a±2, (Xti, 
mi/m2, x\, and <J\j(Ji- For the sake of concreteness, henceforth we will consider the case a\\ = a2i — cti2 = ol. In 
Fig. we show the dependence of 7 on a in the three-dimensional case for Xx = §, o 1 /o"2 = 1 and three different 
values of the mass ratio mi/rri2- We include the analytical results obtained for IMM and IHS2£ as well as the results 
(symbols) obtained^ by numerically solving the Boltzmann equation by means of the direct simulation Monte Carlo 
(DSMC) method^ It is apparent that the analytical results for IMM and IHS are practically indistinguishable over 
the range of values of a considered and that the agreement of both approaches with simulation is excellent. We also 
observe that the extent of the equipartition violation is greater when the mass disparity is large. In addition, the 
temperature of the excess species is larger (smaller) than that of the defect species when the excess species is heavier 
(lighter) than the defect species. The dependence of the coefficients c\ and C2 on the coefficient of restitution a is 
plotted in Fig. Q]in the three-dimensional case for x\ — \, <J\I<J2 = 1 and m\jra2 = 2. It can be observed that the 
HCS of IMM deviates from the Gaussian distribution (which corresponds to c\ = C2 = 0) much more than the HCS 
of IHS. This is consistent with the fact that the former models have a stronger overpopulated high energy tail^iSiiS 
than the latter 4& For both interaction models, the deviation of $ r from its Gaussian form is more significant for the 
heavy species. We also see that the value of the kurtosis c r predicted by the IMM exhibits quantitative discrepancies 
with the one found for IHS, especially for strong dissipation. 



IV. CHAPMAN-ENSKOG SOLUTION OF THE BOLTZMANN EQUATION FOR IMM 

In this Section, the Chapman-Enskog methodic generalized to inelastic collisions will be applied to the set of 
Boltzmann equations Q for IMM to get explicit expressions for the Navier-Stokes transport coefficients as functions 
of the coefficients of restitution and the parameters of the mixture (masses, composition, and sizes). 

The balance equations H12|) - (|14[) become a closed set of hydrodynamic equations for the fields n r , u, and T once 
the fluxes l|15l) - (|17|l and the cooling rate ( are obtained in terms of the hydrodynamic fields and their gradients. 
As noted in Sec. [I] while the pressure tensor has the same form as for a one-component system, there is greater 
freedom in representing the heat and mass fluxes. Here, as done in the IHS casepl we take the gradients of the mole 
fraction x\ = ni/n, the pressure p, the temperature T, and the flow velocity u as the relevant ones. Thus, in this 
representation, the phenomenological constitutive relations for the fluxes in the low-density regime have the formsSS 

mim2fi p p 



ji = 1 -^DVx 1 ~^D p V P ~^D'VT, ja = -ji, (40) 



q = ~T 2 D"Vx 1 - LVp - AVT, (41) 



P t j = pSij - rj yVjUi + \7 iUj - -S i:j ■ V • uj . (42) 

The transport coefficients in these equations are the diffusion coefficient D, the thermal diffusion coefficient D' , the 
pressure diffusion coefficient D pi the Dufour coefficient D" , the thermal conductivity A, the pressure energy coefficient 
L, and the shear viscosity rj. 
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The Chapman-Enskog method assumes the existence of a normal solution in which all space and time dependence 
of the distribution function occurs through a functional dependence on the hydrodynamic fields 



/ r (r,v,t) = f r M Xl (t),p(t),T(t),u(t)}. 



(43) 



This functional dependence can be made local in space and time by means of an expansion in gradients of the fields. 
Thus, we write f r as a series expansion in a formal parameter e measuring the nonuniformity of the system, 



/r=/ f (0) + £ / r (1) + £ 2 /i 2) +--- 



(44) 



where each factor of e means an implicit gradient of a hydrodynamic field. The local reference state fr°^ is chosen to 
give the same first moments as the exact distribution f r , or equivalently, the remainder of the expansion must obey 
the orthogonality conditions 



dv 



/,(v)-/i 0) (v) 



= 0, 



(45) 



J dvm r v [/ r ( v )-/(°)(v) 



= 0, 



(46) 



E / dv ir 2 [Mv)-/i o) (v) 



0. 



(47) 



Use of the expansion Q44|l in the definitions of the fluxes (|15|) - H17|1 and the cooling rate (|10(l gives the corresponding 
expansion for these quantities. The time derivatives of the fields are also expanded as dt — d 



(0) 



ed., 



(i) 



The 



coefficients of the time derivative expansion are identified from the balance equations (|12() - (|14|l after expanding the 
fluxes and the cooling rate £. In particular, the macroscopic balance equations to zeroth order become 



d[ 0) Xr = o, dl 0) u = o, T- l d™T = p-'d^p = -c (0) 



)(<>) 



-la(0)r 



(48) 



Here, we have taken into account that in the Boltzmann equation Q the effective collision frequency ui rs oc n^T 1 / 2 oc 
x s pT~ x l 2 is assumed to be a functional of f r and f s only through the mole fraction x s , the pressure p, and the 



temperature T. As a consequence, w« = w rs , and 



(2) 



0. 



In the zeroth order, fr obeys the kinetic equation 

(V/W)=£J^°U (0) ], (49) 

S 

where use has been made of the relation (|33|1 with = £r given by Eq. The distribution /,^ 0) is given by the 

scaling form l|28l) except that n r — > n r (r, t) and T — > T(r, t) are local quantities and v — > V = v — u(r, t). Since 
is isotropic, it follows that 



I,(o) A 

2^ <9V 



ji 0) = 0, 



q (0) -0, ^ 



(0) 



where p = nT is the hydrostatic pressure. 

In the first order, the distribution function verifies the kinetic equation 



(d^ + Cr) + M r fP = (DP + V • V) /<°>, 



(50) 



(51) 



where it is understood that r ^ s. Here, -Dj 1 " 1 = 9^ + u • V and we have introduced the linearized collision operators 



C 



(52) 



(53) 
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The action of the material time derivatives on the hydrodynamic fields is 



D[ 1] xx = 0, 



d+2 



d[ 1] Inp = -D { t 1] InT = -V • u, L» t (1) u = -p^Vp. 



(54) 



In these equations use has been made of the results (|50f) and the identity £W = 0. The last equality easily follows 

from Eq. I|18|) . The right-hand side of Eq. I|51|) can be evaluated by using Eqs. I|54|l and so the equation for can 
be written as 



(<9 t (0) + £ r ) + MrfP = A r ■ Vxt + B r • Vp + C r • VT + D r>ij V 



(55) 



where 



A r (V) 



(56) 



B r (V) = - 



f(0)y + l ( if(») 

Jr P \av Jr 



(57) 



Cr(V) = - 



V. 



(58) 



D r ,ij(V) 



dVi 



(v,/<»>) - i 



(59) 



It is worth noting that Eqs. I|55 |l -I|59 |l have the same structure as that of the Boltzmann equation for IHS. 21 The only 
difference between both models lies in the explicit form of the linearized operators C r and M. r . 

Now we are in conditions to get the expressions for the mass flux, the pressure tensor, and the heat flux in the first 
order of gradients. These expressions allows one to identify the relevant transport coefficients of the mixture through 
Eqs. |@0H12J}. 

V. NAVIER-STOKES TRANSPORT COEFFICIENTS 

This Section is devoted to the determination of the Navier-Stokes transport coefficients associated with the irre- 
versible fluxes. We only display here the final expressions for the transport coefficients with technical details given in 
Appendix [51 

To first order in the hydrodynamic gradients, the mass flux has the form given by Eq. (|40|) . The transport coefficients 
D, Dp, and D 1 are given by 



D 



m\m,2n \ 2 



p,T 



pc (0) 



{D p + D') 



(60) 



D„ 



m\nT 
pTi 



£ c (o) + L_ 



(61) 



C(o) 

D' = -- — D v , 



(62) 



where 



v = — — /i 2 i(l + ai2) = 3 



dp2 



d n(m\ + m-2) 



h + e 2 



1\02 



1/2 



v a {l + a 12 ). 



(63) 
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Since ji = — j 2 and Vxi = — Va; 2 , D should be symmetric while D p and D 1 should be antisymmetric with respect to 
the exchange 1 <-> 2 . This can be easily verified by noting that n{T\ + n-iTi — nT. For elastic collisions, = 0, 
Ti=T 2 = T, and so Eqs. ffi ty -ffT ^ become 

D= dm 1+m2 p_ t 
8 m\m2 Pa 



D P = a K m 2 ~ m i)— 3 £>=0, 

8 p- l/Q 



(65) 



where use has been made again of Eq. 1)22)1 . In the case of mechanically equivalent particles, 7=1, D p = D' = 0, so 
that Eq. JHDJ gives the expression of the self-diffusion coefficient 



m \ 2 

The pressure tensor has the form i|42[) with the shear viscosity coefficient rj given by 

r/ = i]i+ i] 2 . 

The expression of the partial contributions r\ r is 

Pi(2r 22 -C (0) )-2p 2 r 12 



m = 2 



C (0)2 _ 2C(°)(r 11 + r 22 ) + 4(ni72i 2 - r 12 r 21 ) : 



(66) 



(67) 



(68) 



where the quantities m and n 2 are defined by Eqs. I)B12)I and IB 1 3)1 . respectively. A similar expression can be obtained 
for rj2 by just making the changes 1 «-> 2. For mechanically equivalent particles, Eq. (|68|l yields r\\jx\ = 77 2 /x 2 = r/, 
where 



n 



p 



- I((o) 



(l + a)(rf+ I -a) 
d{d + 2) 



-uj, 



(69) 



and uj = uj rs /x s . The expression (IrJ^t coincides with the one previously derived in the single gas case^ 

The case of the heat flux is more involved. Its form is given by Eq. I|41|l where the coefficients D", L and A are 



D" = D'{ + D'l L = Li + L 2 , A = Ai + A 2 . 
By using matrix notation, the coupled set of six equations for the unknowns 

{D 1 , D 2 ,Li,L2, Ai, A 2 } 

can be written as 

A-aa' 'X a f Y a . 

Here, X a i is the column matrix defined by the set (|71|l and K aa i is the square matrix 



(70) 
(71) 
(72) 







A 



V 



-T 2 /3 21 








T 2 (§C (0) - 022) 







T 



\ 9x1 ) P a 



\ 9x1 J P ,i 







T 



-Pi 



22 



The column matrix Y is 



C (0) -/3n 

-021 

-K (0) /2T 

-p(^/2T 



Y 2 
Y 3 
Y 4 
Y 5 

\yJ 



T((°)/p 


C (0) - Pa 

-021 



/ ac (0) 





-012 
C iO) - 022 J 



(73) 



Y 



(74) 
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where 



Tn\Tn 2 n . d + 2nT 2 d 

Y\ = A12-D 7- — 

p 2 mi t> x i 



1 + |) 



(75) 



rti\m 2 n d + 2 nT 2 d 

Y 2 = A 2 \D — — 

p 2 m 2 0x1 



1 + f ) 



(76) 



Y 3 = ^A 12 D P 



d + 2 n{Tl / _ mip £1 
2 mip \ pTi 2 J 



(77) 



d + 2 n 2 T 2 ( \ m 2 p , c 2 



Y 4 = -^A 21 D p 

p 2 m 2 p \ pl 2 2 



(78) 



Y 5 = ^A 12 D 



, d + 2 m2f 
2 miT 



Cl 



(79) 



The expressions of the quantities f3 r 
The solution to Eq. is 



+ 2 n 2 T 2 
~2~ m 2 T 



! + -=• 



and A rs are given in Appendix IbI 



X a = {K-') aal Y a 



(80) 



(81) 



This relation provides an explicit expression for the coefficients D", L r and A r in terms of the coefficients of restitution 
and the parameters of the mixture. From these expressions one easily gets the transport coefficients D" , L, and A 
from Eq. I|7U|) . As expected, Eqs. (|72|I - H8U|I show that D" is antisymmetric with respect to the change 1 <-» 2, while 
1/ and A are symmetric. This implies that in the case of mechanically equivalent particles, the Dufour coefficient D" 
vanishes. Furthermore, in this limit, Eq. (|81|l leads to the following expression for the heat flux 



q= -kVT - pWn, 



(82) 



where 



k = A + nL = 



d + 2 p 1 + c 
2 mv K - 2CW ' 



(83) 



/i = TL = 



T k(C (0) + 5C^) 
n(l + c )K-|C( 0) ) : 



(84) 



with 



4(d-l) + (8 + d)(l-a) 
4d + 4(1 - a) ' 



(85) 



Note that upon writing Eq. I|82|l use has been made of the relation Vp = nVT + TVn. Again, Eqs. (|82() - (|85|l coincide 
with results derived for a single gas4& Using Eqs. (|83|l Q85J1 ■ it can be seen that the coefficients k and p diverge for 
d — 2 and d = 3 when a = ao = (4 — d)/3d and so, both coefficients become negative for < a < olq. This unphysical 
behavior could be due to the failure of the hydrodynamic description for values of the coefficient of restitution smaller 
than ao or perhaps to the existence of a certain type of hydrodynamic instability. However, given that the value of ckq 
is quite small (ao — -| at d — 2 and cto = ^ at d — 3) this singular behavior could be interpreted as a hydrodynamic 
breakdown indeed. Elucidation of this point requires further analysis. It must be remarked that the above drawback 
is absent for IHS since all the transport coefficients are regular functions of a for all d. 
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a 



FIG. 3: Plot of the reduced self-diffusion coefficient D(a) /D(l) as a function of the coefficient of restitution a in the three- 
dimensional case for IMM (solid line) and for IHS in the first Sonine approximation (dashed line)— and in the second Sonine 
approximation (dotted line)i— ■ The symbols refer to Monte Carlo simulations for IHSi— *— 




0.5 0.6 0.7 0.8 0.9 1.0 

a 



FIG. 4: Plot of the reduced diffusion coefficient D(a)/ D(l) as a function of the coefficient of restitution a in the three- 
dimensional case for <7i/<72 = 2, x\ — (tracer limit) and mi/7712 = 8. The solid line corresponds to the exact results obtained 
here for IMM while the dashed line is the result derived for IHS in the second Sonine approximation^ The symbols refer to 
Monte Carlo simulations for IHSp— 

VI. COMPARISON WITH THE TRANSPORT COEFFICIENTS FOR IHS 

The expressions for the transport coefficients of IHS described by the Boltzmann equation have been obtained 
by Garzo and Dufty^i in the leading Sonine approximation for a three-dimensional system. These expressions have 
been then evaluated for a variety of mass and diameter ratios in the cases of the diffusion coefficient and the shear 
viscosity coefficient, 2 - showing quite a good agreement with Monte Carlo simulations. In this section, we compare 
the results derived here for IMM for the transport coefficients entering in the mass and momentum fluxes, namely, 




a 



FIG. 5: Plot of the reduced diffusion coefficient D(a)/ D(l) as a function of the coefficient of restitution a in the three- 
dimensional case for ci = <72, xi = 0.2, and two different values of the mass ratio: mi/7712 = 0.5 and m\/m2 = 4. The solid 
lines correspond to the exact results obtained here for IMM while the dashed lines are the results derived for IHS in the first 
Sonine approximationi— 
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a 



FIG. 6: Plot of the reduced pressure diffusion coefficient D p (a) / D p (l) as a function of the coefficient of restitution a in the 
three-dimensional case for a\ = a%, x\ = 0.2, and two different values of the mass ratio: mi/ma = 0.5 and m\jm-i — 4. The 
solid lines correspond to the exact results obtained here for fMM while the dashed lines are the results derived for IHS in the 
first Sonine approximation!— 




a 



FIG. 7: Plot of the reduced coefficient —D*(a) = — D' (a) / D p (l) as a function of the coefficient of restitution a in the three- 
dimensional case for oi = <72, x\ = 0.2, and two different values of the mass ratio: m\/m2 = 0.5 and mi/m2 = 4. The solid 
lines correspond to the exact results obtained here for fMM while the dashed lines are the results derived for IHS in the first 
Sonine approximation!— 

D, Dp, D', and r\ with those obtained for IHS in the case d — 3. As in the cases studied in Section ITTT1 we consider for 
the sake of simplicity a common coefficient of restitution, i.e., an = a\2 — (X22 = a. 

Let us consider first, the diffusion coefficient D in the limit cases of self-diffusion (mechanically equivalent particles) 
and tracer concentration (x± — > 0). Figure |21 shows the reduced self-diffusion coefficient D(a)/ D(l) as a function 
of the coefficient of restitution a. Here, D(l) refers to the self-diffusion coefficient for elastic collisions. We 
include the results obtained for IHS by using the first Sonine approximation (dashed Iine)j2i2i the second Sonine 
approximation (dotted line) and by Monte Carlo simulations (symbols).— We observe that the agreement between 
the predictions of the first Sonine approximation for IHS and for IMM is excellent in the whole range of values of 
a analyzed. Moreover, both theories compare quite well with computer simulations even beyond the quasielastic 
limit (say for instance, a > 0.8). However, as dissipation increases, the agreement between theory and simulation 
is improved when one considers the second Sonine approximation. For mechanically different particles, in Fig. 0] we 
plot D(a) I D(l) versus a in the tracer limit (xi — > 0) for mi/m,2 — 8 and (T\ja% = 2 as given by the exact results 
for IMM, the second Sonine approximation for IHS 26 and by Monte Carlo simulations. We see that both theories are 
practically indistinguishable and present an excellent agreement with simulation data. 

Beyond the above two special cases, Figs. show the dependence of the reduced coefficients D(a)/D(l), 
D p (a)/Dp(l), and — D*(a) on the coefficient of restitution for a\ — 02, x\ — 0.2, and two values of the mass ratio 
(mi/ma = 0.5 and 4). Here, D p (l) is given by the first equality of Eq. (HSl and D*(a) = -(C (0) /2v)D p {a) / D p (l) = 
—D'(a)/D p (\). The solid lines correspond to the IMM results while the dashed lines refer to the results obtained for 
IHS in the first Sonine approximation^! We observe that, in general, the qualitative behavior of IHS is well captured 
by the IMM. At a quantitative level, for not strong dissipation (say for instance, a > 0.8) the agreement between 
the results derived for IMM and IHS is again quite good, especially in the case of mi/m<z = 0.5. Nevertheless, the 
discrepancies between both interaction models increase as the coefficient of restitution decreases. For instance, at 
a = 0.5, the discrepancies for D(a)/D(l), D p (a)/D p (l), and — D*(a) are about 4%, 5%, and 8%, respectively, for 
TOi/m 2 = 4, while they are about 1%, 3%, and 3%, respectively, for mi/m 2 = 0.5. 
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FIG. 8: Plot of the reduced shear viscosity coefficient r)(a)/r](l) as a function of the coefficient of restitution a in the three- 
dimensional case for <7i = a%, xi = 0.2, and two different values of the mass ratio: m\/rri2 = 0.5 and mi/mz = 4. The solid 
lines correspond to the exact results obtained here for IMM while the dashed lines are the results derived for IHS in the first 
Sonine approximation!— 

The dependence of the reduced shear viscosity coefficient rj(a)/r](l) on dissipation is plotted in Fig. for a\ = <t 2 . 
x\ = 0.2, and for two values of the mass ratio {mi/mi = 0.5 and 4). Here, 77(1) is the corresponding value of the 
shear viscosity coefficient for elastic collisions. Although again the IMM predictions compare qualitatively well with 
the ones derived for IHS, we observe that in general the discrepancies for both interaction models at the level of the 
shear viscosity are larger than those found for the transport coefficients of the mass flux. The discrepancies observed 
in Fig. |8] are similar to those found in the single gas case^ 

VII. DISCUSSION 

The primary objective of this work has been to derive the hydrodynamic equations of a granular binary mixture from 
the Boltzmann kinetic theory for inelastic Maxwell models (IMM). In the Boltzmann equation for inelastic Maxwell 
models (IMM), the collision rate of inelastic hard spheres (IHS) is replaced by an effective collision rate independent 
of the relative velocity of the two colliding particles. As in the elastic case^S this property allows us to get exactly 
the velocity moments of the Boltzmann collision integrals without the explicit knowledge of the velocity distribution 
function. Here, the Chapman-Enskog method has been applied to get a normal solution of the Boltzmann equation 
for states near the local homogenous cooling state. The derivation of the Navier-Stokes hydrodynamic equations 
consists of two steps. As a first step, the reference homogenous cooling state for a mixture of inelastic Maxwell gases 
is analyzed to provide the proper basis for description of transport due to spatial inhomogeneities. As in the case 
of inelastic hard spheres (IHS)j22 our solution for the homogeneous state shows that the kinetic temperatures for 
each species are clearly different so that the total energy is not equally distributed between both species (breakdown 
of energy equipartition). In addition, we also compute the fourth cumulant of the velocity distribution functions, 
which is a measure of the deviation of the distributions from the Maxwellian form. Once the reference state is well 
characterized, as a second step we obtain exact expressions for the mass flux, the pressure tensor, and the heat flux 
in the first order of the hydrodynamic gradients (Navier-Stokes order). From these expressions we identify the seven 
relevant transport coefficients of the problem, namely the mutual diffusion D, Eq. I|60[) . the pressure diffusion D p , 
Eq. (|ulf> . the thermal diffusion D', Eq. (|62|l . the shear viscosity 77, Eqs. I|67l) and (|68[) . the Dufour coefficient D", 
the pressure energy coefficient L, and the thermal conductivity A as given by Eqs. 170|) l)8ip. respectively. These 
expressions are exact (within the context of IMM) and constitute the main goal of this paper. This contrasts with the 
previous results derived for IHS, 21 where the transport coefficients have been approximately obtained by considering 
the leading terms in a Sonine polynomial expansion of the distribution function. 

The purpose of our work is twofold. First, the evaluation of the transport coefficients for mixtures of IMM is worth- 
while studying by itself as a simple model to gain some insight into the influence of collisional cooling on the transport 
properties of the system. Second, the comparison between the exact results for IMM with the results obtained for 
IHS by using approximate analytical methods allows us to assess the degree of reliability of IMM to reproduce the 
relevant behavior of IHS in the context of granular mixtures. Recent results^ derived for multicomponent systems in 
the simple shear flow problem show a good agreement between both interaction models over a wide range values of 
the parameter space. 

To make contact with the results obtained for IHS^i one needs to fix the collision frequencies ui rs . These quantities 
can be seen as free parameters of the model to be adjusted to optimize the agreement with IHS. As in the study made 
in Ref. HH , here we have chosen ui rs to reproduce the cooling rates <^ rs of IHS in the local equilibrium approximation, 
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FIG. 9: Plot of the reduced coefficients «(a)//t(l) and p* = np/T«(l) as a function of the coefficient of restitution a in the 
three-dimensional case for a single gas. The solid lines correspond to the exact results obtained for IMM— while the dashed 
lines are the results derived for IHS in the first Sonine approximation!— 



Eq. i|21|) . An exploration of the dependence of the seven transport coefficients on the full parameter space (mass 
ratio, diameters, concentrations, coefficients of restitution) is straightforward but perhaps beyond the scope of this 
presentation. Here, we have focused on the dependence of the coefficients D, D p , D' , and 7/ on the common restitution 
coefficient a for different values of composition, sizes, and masses. The comparison with known results for iHSS2*2SiS& 
has been illustrated in Figs. |3HH| In addition, Monte Carlo simulations have been also included in some plots. 
The comparison shows that, in general, the IMM predictions are reasonably good for not too large dissipation (say, 
a > 0.8), especially for the transport coefficients associated with the mass flux (first-degree velocity moment). The 
discrepancies between both interaction models increase in the case of the shear viscosity coefficient, which is related to 
a second-degree moment (pressure tensor). However, the IMM results still capture qualitatively well the dependence 
of r] on a for IHS since the discrepancies between both models are about 4% at a = 0.8 and about 9% at a — 0.5. As 
in the monocomponent case^& more significant disagreement between IMM and IHS are expected when one compares 
higher degree moments, such as the heat flux. To show these discrepancies, in Fig. we plot the reduced coefficients 
k(q)/k{\) and p* = np/Tn(l) for d = 3 in the single gas case. Here, k(o) is given by Eq. I|83|l while the coefficient 
/1 (a) is given by Eq. (|84|) . As was emphasized in Ref. flfi the trends observed for IHS are strongly exaggerated by the 
IMM, where k and p diverge at a = |. 

A simple application of the results obtained in this paper would be the determination of the dispersion relations 
for the hydrodynamic equations linearized about a homogeneous state. This analysis would allow us to identify the 
conditions for stability as functions of the wave vector, the dissipation, and the parameters of the mixture. Another 
possible direction of study would be to check the validity of Onsager's reciprocal relations 20 " between the different 
transport coefficients associated with the mass and heat fluxes. Since the system is not time reversal invariant, 
violation of Onsager's relations is expected for inelastic collisions. The goal would be to assess the influence of 
dissipation on such violation. We plan to work on these issues in the next future. 



APPENDIX A: COLLISIONAL MOMENTS IN THE INELASTIC MAXWELL MODEL 

In this Appendix we will determine the collisional integrals of m r V, rrirVV, and m r V 2 'V appearing in the evaluation 
of the transport coefficients. The two first integrals were already obtained by one of the authors in a previous papen^ 
on the analysis of rheological properties. Now, for the sake of completeness, we display their explicit expressions: 



dvm r VJ rs [/ r , / s ] = -/J sr (l + a rs ) (p s j r - p r j s ) , 

p s d 



(Al) 



dvm r WJ rs [f r , f a 



Psd 
2 



+ 



d + 2 
d 



p sr {l + a rs ) {2p s P r - (j,.j s + j s j r ) 
p, sr (l + a rs ) [p s P r + p r P s - (j r j s + j s j r ) 



(p r p s + PsPr) ~jr h 



(A2) 
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where 1 is the d x d unit tensor. It only remains here to get the collisional integral Q rs corresponding to the heat 
flux: 

Q JS = J dv^m r V 2 VJ rs [v\f r ,f s }. (A3) 
To simplify its calculation, a useful identity for an arbitrary function h(y) is given by 

J dvifc(vi)J P ,[v 1 |/ r ,/ a ] = ^ / dvi j dva/r(vi)/,(v 2 ) j da [h(v'{) - h(yx)] , (A4) 

with 

v'/ = V! - fi sr (l + a rs )(a ■ gi 2 )a, (A5) 
and gi2 = Vi — v 2 . Now, we particularize to h(V) = \m r NV 2 . In this case, using 1A4() one has 

Q, s - J dvi | dv a / r (v!)/ a (v a ) | dS (vpV'{ If Vt) . (A6) 

From the scattering rule i|A5[l it follows that 

a V? - V?Va = -fx sr (l + a rs )(a • g 12 ) { [V? - 2 Msr (l + a rs )(a ■ g 12 )(S • V x ) 
+^ r (l + a rs ) 2 (5 ; -g 12 ) 2 ] 

+ [2(5 • Vi) - M S r(l + a rs )(a • g 12 )] Vr} . (A7) 
To perform the angular integration, one needs the results 

J da {a- gi 2 ) fc ? = B k+ ig\z l % 12 , (A8) 
j da {a ■ % 12 ) k aa = ^-/^ (fcgi2gi2 + g\ 2 t) , (A9) 



where 4 



,._„^r(j)rm 



b„ = jd»{9- g 12 )« = n^r-" r ( M) ■ (A10) 

Taking into account Eqs. I|A7(I and l)A8|l - l|A10jl . the integration over 5 in Eq. (|A6(I leads to 

Qrs = ~ - ^ ^Y^sr{l + a rs ) y dvi y dv 2 / r (vi)/ s (v 2 ) 

x [(d + 2)Vfg 12 - 4/x sr (l + a rs ) (gi 2 ■ Vi) gi 2 - (d+ 4)// sr (l + a rs ) 5 2 2 Vi 

+3/4(1 + a rs ) 2 3 2 2 g 12 + 2(d + 2) (g 12 • Vi) Vi] . (All) 

The corresponding integrations over velocity give the relations 

dwx j dv 2 — Vt gi2/r(vi)/ a (v 2 ) = n s q r - T^Ms, ( A12 ) 

y dvi y dv 2 ^(gi 2 • Vi)gi 2 / r (vi)/ s (v 2 ) = n s q r - — ^- (dp r j s + P r ■ j s - P s ■ j r ) , (A13) 

y^vi y dv 2 ^g 2 2 Vi/ r (vi)/ s (v 2 ) = n s q r + — ^- (dp s j r - 2P r • j s ) , (A14) 
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J dvi J rfv 2 ^(gi2 ■ Vi)Vi/ r (vi)/ a (v 2 ) = n s q r - T^-FV js, 

/r jjj xrt 
rfvi / dv 2 — -g 2 2 gi 2 /r(vi)/ s (v 2 ) = n s q r -n r qr 

- ^ (<Zp r j a - rfPsjr + 2P r • j s - 2P S • j r ) 



(A15) 



(A16) 



From Eqs. l|A"T2")) -l|A"T5 )l one finally gets the expression of Q, 



p s d(d + 2) 



(1 + a sr ) {[/t sr (l + a rs ) (<2+ 8 - 3/t ar (l + a rs )) - 3(d + 2)] p s q r 



+3^ r (l + a sr ) 2 p r q s + - \p sr (l + a rs ) (3/t sr (l + a rs ) - 4) + d + 2}p r j s 

(1 + a rs ) [d + 4 - 3/x sr (l + a rg )] p s j r 

+ [m<w(1 + "rs) (3jU sr (l + a rs ) - (d + 6)) + d + 2} P r ■ j s 
(1 + a rs ) [2 - 3/t sr (l + a rs )] P s • j r } , 



(A17) 



where p r — n r T r . In the absence of diffusion and for mechanically equivalent particles, the collisional moment Q rs = Q 
reduces to 



{d- 1) , 



1 - a 



d-1 4 



(A18) 



This expression coincides with the one previously derived in the monocomponent case^& 
Finally, let us evaluate the (dimensionless) collision integral of v 4 in the HCS, Eq. I|37|) : 



A rs = J dv*«* 4 J r * s [$ r ,$ a ] 

dv^ [ da<$> r {vl)<$> s (v* 2 ) (v'{ 



l/i 1 



(A19) 



where use has been made of the property (|A4J| . Henceforth, we will use dimensionless quantities and, for the sake of 
simplicity, the asterisks will be deleted. The scattering rule (|A5|) gives 



v'f-v\ = 2/4(1 + a rs f {* -g 12 ) 2 2(?. Vl ) 2 + 2^ + /|n(i + aj , s )2 ( ^.g 12 )2 
-4/i ar (l + • g X2 )(CT ■ vi) [u 2 + p 2 sr (l + a rs f{a ■ gi 2 ) 2 ] . 

Equations lA8j) - l)A10jl allows the angular integral to be performed with the result 



J dcT^-vf) 



4B 
~d+2 

4B 2 



2( Vl • g 12 ) 2 + ^^g 2 12 v 2 + ^ 2 r (l + a rs ) 2 gi 2 



2 p 2 sr {l + a rs ) 2 

/t sr (l + a rs )(vi • gi 2 ) [(d + 2)v\ + 3// 2 r (l + 



d+2 



Therefore, the collision integral A rs can be written as 



A rs = -^-B 2 p 2 sr (l + a rs ) 2 I dvi I dv 2 $ r (ui)$ s (w 2 ) 



1 



+3 



d + 2 
p 2 r (l + a rs ) 2 
d + 2 



2{d + 8) - 12/i sr (l + a rs ) + 3/4(1 + a rs ) 2 - U " + 2 ' 



fi sr (l + a rs ) 



2 

v i + -[d + 2- 6/x sr (l + a rs ) + 3/4(1 + a rs ) 2 ] v 2 v 2 



(A20) 



(A21) 



(A22) 
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Finally, by taking into account that 



/d 



dvv^ r (v) = ( A23 ) 



we get 



[i sr (l + a rs ) , 3 3 4 



+ [2d + 3^(1 + a rs ) 2 - 6/i sr (l + a rs ) + 4] [/i sr (l + cv s ) - 2] («* 4 ) r 
-/^"^(l + a rs ) [2d + 4 + 6^ sr (l + a rs ) 



4 

x (^ sr (l + a rs ) -2)]}. (A24) 

In the one-dimensional case (d = 1), Eq. (|A24|I agrees with the results derived in Refill for the scalar Maxwell model 
(i.e., with uj* s oc x s ). In addition, for mechanically equivalent particles, the results of the single gas are recovered, 
namely^ 

A rs = A = 23d( ^ + (1 + a) { [I2a 2 (a - 1) + 4a(4d + 17) - 12(3 + 4d)( v 4 }] 

+ d(d + 2) (i + Q)(4rf _ 16a + 3q2) 1 _ (A25) 



APPENDIX B: DERIVATION OF THE TRANSPORT COEFFICIENTS 

In this Appendix we will provide some details on the calculation of the Navier-Stokes transport coefficients appearing 
in the expressions ljin jl -l(4"2* jl of the irreversible fluxes. Let us consider each flux separately. 

1. Mass flux 



To first order, the mass flux j[ is defined as 



i^ = mi fdvyfP(v). (bi) 

To get this flux from Eq. (|55|l . we need the collisional integral of m{V which has been evaluated in the Appendix lAl 
From the linearization of Eq. (|A1|) , one has the result 

J dv mi V (£ifl X) +M 1 f£ ) ) = v'£\ (B2) 

where v is the collision frequency 

no. ^ 4 p ( 6 l+ 6 2 \ 1/2 
dp 2 dn(mi+m 2 j \ 0i0 2 J 

and use has been made of Eq. lt%2l in the second equality. Next, we multiply both sides of Eq. by miV and 
integrate over V. The result is 

d \ Ti\Ti ( ^ mip 



KHtf' = -> stH. t '-T '-i^ (B4) 



Note that the temperature ratio 71 depends on the hydrodynamic state through the concentration X\. The functional 
dependence of 71 on X\ can be obtained from the HCS condition l|29|l by using the expressions (|27H for the partial 

cooling rates Cr°' ■ 
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The mass flux has the structure given by Eq. Q40[l. Dimensional analysis shows that D oc T 1 / 2 , D p oc T^/p, and 
D' oc T 1 / 2 . Consequently, 



p V 2 r ' ) r 2T p 
Upon deriving this expression use has been made of the identities 



(0) 



VT = -V (TC (0) ) = -C (0) VT-TVC (0) 



£(0) 



VT-T 



/0 C (o)\ C (o) 



(B5) 



(B6) 



a t (0) vp = -v(pc (0) ) =-c (0) vp-pVC 



where we have taken into account that 



(0) 



2C (0) Vp - p 



dx i 



f(0) 



V / p T V Op / xi, T 

/ 5C (o)x C (o) C (0) 



G>C (0) 

dT 



VT 



(B7) 



(B8) 



the two latter terms coming from Q ' oc nT 1 ! 2 = pT 1 / 2 . Inserting Eq. (|B5J| into Eq. I|B4(I . one gets the expressions 
1|60J) . (|61|) . and (|62|) for the coefficients D, D p , and D', respectively. 



2. Pressure tensor 



The pressure tensor pW can be written as 



p(i) = P W +p (i) ) 



where the partial contribution Pr to the pressure tensor is 

PW =m r J dvVV/W(v). 
The linearization of Eq. (|A2|) leads to the following expression for the collisional integral of miVV: 



J dvrnxVV (A/W +MifP) = ruPW +r 12 pW, 



where 



^11 ^12 

T n = , I 1 + + 1 - + 2— ^2l(l + 012) 



d(d + 2) 



1 - 



M2l(l + «12) 

d + 2 



^12 Pi 2/i 1 \2 

Tl2 = ~ 2 Tv77T^ — ^21 1 + "12) • 
d(d + 2) p 2 



(B9) 
(BIO) 

(Bll) 

(B12) 
(B13) 
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Now, we multiply both sides of Eq. i|55|l (with r = 1) by miVV and integrate over V to get 

(d t (0) + r n ) P[% + T 12 P^l = - Pl A t3ke V k u e , (B14) 



where p\ — n\T\ and 



2 

A-ijM = SikSje + SitSjk — —SijSke- (B15) 



A similar equation can be obtained for P 2 from l)B14jl by interchanging 1^2. The solution to Eq. 1B14J1 (and its 
corresponding counterpart) has the form 

P$j = -Vr&ijktVkUi. (B16) 
According to Eq. I|42(l . the shear viscosity coefficient r\ is given in terms of the coefficients r\ r by 

77 = 771+772- (B17) 
Dimensional analysis requires that ?7 r oc T 1 / 2 and so, 

9 (o)p(i) = _C^p(D_ (Blg) 

Insertion of this relation into Eq. I|B14|) yields the following set of coupled equations for the two coefficients rj r : 

ru-|C (0) r 12 _,\fvi\ = f P i\ (B19) 



T21 T 22 - iC (0) J \ V2 J \P2 

Its solution is given by Eq. . 

3. Heat flux 

The heat flux can be written as 

q^qW+qW (B20) 



where the partial contribution q£^ is given by 



(1) m ! 
^ 2" 



/ 'dvV*VfP(v). (B2I) 



To get the explicit expressions for the fluxes q^ we proceed in a similar way as in the case of the pressure tensor. 
First, linearization of Eq. (|A17I) leads to 

/ d,^V (CJ^+MJ^) = AxqW +/3 12 q 2 1) + A 12 j« (B22) 



where 



Wn(l + o; u ) r ,,. e s ,1 (l + ai2) 

x {^21(1 + a u ) [d + 8- 3M2i(l + o 12 )] - 3(d + 2)} , (B23) 



ft - % , ..3 (l + «12) 3 Pl 

P12 - -Aji 2 /x 21 , (B24) 

d(c! + 2) p 2 
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T 2 



x ^2i(l + ai 2 )[d-3At2i(l + ai 2 ) + 2] — 



— [d + 3/4(1 + a 12 ) 2 - 6 M 2i(l + au) + 2] — L 
z 2 m 2 1 



(B25) 



Upon writing Eq. I|B22|1 use has been made of the relation = — j| . The corresponding expressions for /3 22 , /3 2 i, 
and A21 can be easily obtained from Eqs. I|B23(I l|B25() by the change 1 <-> 2. From Eq. I|55|) . one gets 



(^° ) + /3 11 )qi 1) +A 2 q 2 1) = -A 12 j«- 



d + 2nT 2 d 



2 mi <9xi 



d + 2 mTf 

2 TOxp 

d + 2 mTf 
2~ miT 



1 



pTi 2 
x VT, 



1+f )• 



where 



d(d + 2) n r T 2 



The coefficients c r have been obtained in Section ITTT1 
The solution to Eq. (|B26|) can be written as 



q« = -T 2 D'y Xl ~ L r Vp - A r VT. 



(B26) 



(B27) 



(B28) 



The total heat flux defines the transport coefficients D", L, and A through Eq. I|4ip. According to Eqs. I|B20|) and 
(|B28|) . these transport coefficients are given in terms of their partial contributions D", L r) and X r as 

D" = D'l + D'2, L = L X + L 2 , A = Ai+A 2 . (B29) 

From dimensional analysis, D" cx T -1 / 2 , L r cx T 3 / 2 /p, and \ r cx T 1 / 2 . Consequently, 

3 



9 t (0) q« = 



+C (0) (l^r + ^Vp + C^ (A, 



2T 



Van 
VT. 



(B30) 



Substitution of Eq. i|B30(l into Eq. (|B26|I and taking into account the expression Q40f) for the mass flux, one arrives at 
the coupled set of equations H72[) for the partial contributions D", L r , and A r . 
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